Field-induced decay dynamics in square-lattice antiferromagnet 
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Dynamical properties of tlie square-lattice Heisenberg antiferromagnet in applied magnetic field 
are studied for arbitrary value S of the spin. Above the threshold field for two-particle decays, the 
standard spin-wave theory yields singular corrections to the excitation spectrum with logarithmic 
divergences for certain momenta. We develop a self-consistent approximation applicable for S > 
1, which avoids such singularities and provides regularized magnon decay rates. Results for the 
dynamical structure factor obtained in this approach are presented for S — 1 and S = 5/2. 
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I. INTRODUCTION 

The square-lattice Heisenberg antiferromagnet is an 
important model system in quantum magnetismii Theo- 
retical studies of this model have provided deep insights 
into the role of low-dimensionality in the static and dy- 
namical properties of the many-body spin systems. Re- 
cently, there has been a growing interest in the effect 
of magnetic field in the dynamics of quantum antifcrro- 
magnets. The strong-field regime is now reachable for a 
number of layered square-lattice materials with moderate 
strength of exchange coupling between spins3i^— In addi- 
tion, new field-induced dynamical effects can be present 
in the antiferromagnets with other lattice geometries^ 
as well as in the gapped quantum spin systems near the 
condensation field for triplet excitations 

The ground-state properties of the square-lattice anti- 
ferromagnet in a finite field conform with the semiclas- 
sical picture of spins gradually tilted towards the field 
directiouji^ sec Fig. [T] On the other hand, excitation 
spectrum and dynamical properties are expected to un- 
dergo rather dramatic transformation^^"— Following the 
prediction of the field-induced spontaneous magnon de- 
cays in Ref. [H, there is an ongoing search for suitable 
spin- 1/2 square- lattice materials^"— to investigate such 
an effect. The existence of substantial damping in the 
magnon spectrum of the square-lattice Heisenberg model 
in a field has been recently verified by the Quantum 
Monte Carloii (QMC) and the Exact Diagonahzationi^ 
numerical study. Other theoretical aspects of the be- 
havior of the quantum square-lattice antiferromagnet 
(SAFM) in applied field have also been addressed li^Ti^i 

In this paper, we extend previous work of two of usi^ 
and provide a comprehensive theoretical investigation of 
the dynamics of the nearest-neighbor Heisenberg SAFM 
including detailed calculation of the 1 /S correction to the 
energy spectrum in external field, kinematic analysis of 
the field-induced two-magnon decays, and self-consistent 
treatment of magnon decay rates for systems with S > 1. 




FIG. 1: (Color online) Canted spin structure of the square- 
lattice Heisenberg antiferromagnet in external magnetic field. 



The model spin Hamiltonian is 



(1) 



where J stands for the nearest-neighbor exchange cou- 
pling constant and H is the external magnetic field in 
units of g/iB . 

The standard spin-wave theory works quite well for 
the SAFM in zero iieldJ^^— Surprisingly enough, in 
high magnetic fields one encounters strong singularities 
in the spin-wave corrections to the dynamical properties 
of the SAFM, which arise due to spontaneous magnon 
decays above the threshold field H* w 0.75HsJ^ Singu- 
lar behavior of the spectrum signifies a breakdown of the 
perturbative l/S expansion and requires a regularization. 
The aim of this work is to develop a self-consistent ap- 
proximation which yields the spectrum that is free from 
the essential singularities. 

Within the conventional spin-wave approach, the gap- 
less Goldstone branch is preserved only if all quantum 
corrections to the spectrum of the same order in l/S 
are taken into account. This represents an obvious diffi- 
culty for any self-consistent calculation, which typically 
involves summation of a certain infinite subset of pertur- 
bation series. The basic idea of the present work is to ne- 
glect the real part of the spectrum corrections completely 
and to perform the self-consistent calculation only in the 
imaginary part of the magnon self-energy. Such an ap- 



proach is justified if the real part of corrections is small, 
which is the case for S* > 1 for the model ([1]). Utiliz- 
ing this self-consistent scheme we obtain explicit results 
for the magnon decay rates and the dynamical structure 
factor of SAFMs with S = 1 and S = 5/2. 

The paper is organized as follows. Technical details 
of the spin- wave theory in magnetic field are provided in 
Sec. II together with explicit calculation of the quantum 
correction for the spin-1/2 SAFM. The kinematic anal- 
ysis of the field-induced magnon decays and associated 
singularities is presented in Sec. III. The self-consistent 
theory is described in Sec. IV, where the results for the 
dynamical structure factor are also included. Sec. V con- 
tains our conclusions. The asymptotic form for the decay 
rate of low-energy magnons is derived in Appendix A and 
an extension of the kinematic analysis to anisotropic an- 
tifcrromagnets is given in Appendix B. 



II. SPIN- WAVE EXPANSION 

In this section we provide necessary details of the stan- 
dard spin-wave formalism as applied to two-sublattice 
antiferromagnets in external magnetic field at zero 
temperature fi^ii^iSi The first essential step is to quantize 
spin components in the rotating frame {xi,yi,Zi) such 
that the local axis points in the direction of each mag- 
netic sublattice. In the case of the square-lattice model 
the ground state in zero magnetic field is the collinear 
Neel order characterized by the wave- vector Q = {tt,tt). 
In a finite magnetic field spins cant towards the field di- 
rection by an angle 9, see Fig.[TJ Spin components in the 
laboratory frame {xo,yoT zq) are related to those in the 
local frame by 



$1 



5f ° = 5f sin 9 + 5f e''^""' cos 9 , Sf" 
5'f ° = -Sfe'^-''^ cos 9 + Sf sin 9 . 



Sf 



(2) 



The above representation allows us to introduce only 
one type of bosons via the standard Holstein-Primakoff 
transformation<21 Substituting ^ into Eq. ([T]) and ex- 
panding square-roots one obtains bosonic Hamiltonian 
as a sum H = J2'^=o ^n, each term being of the n-th or- 
der in bosonic operators and carrying an explicit factor 
52-n/2^ For large spins this form provides the basis for 
the 1/5* expansion. 

Minimization of the classical energy Hq determines the 
canting angle 



sin 61 = H/Hs with = 8JS 



(3) 



The linear terms given by "Hi vanish for this choice of 9 
and expansion begins with the quadratic terms 

7^2 = H sm9^ala, + JS^^cos29 {ala, + ajaj) 

(ij) 

+ sm^9 {a\aj + aja;) - cos^6l (a|aj -t- ajai) . (4) 




-ft 



k >^ 

k-q+C 




Q-k-g 



FIG. 2: (Color online) The lowest-order normal self-energies 
generated by the decay (left) and by the source (right) vertex. 



Performing successive Fourier and Bogolyubov transfor- 
mations, the latter being defined as 

we obtain the diagonal form of 'H2'- 

?i2 = 5]ek5[5k + i5](6k-Ak) (5) 



with 



ek = 4J5v/(l+7k)(l- cos 207k) = 4J5wk , 

2/1 \ 2 2 ± fk 



(6) 



Ak = 4J5'(l + sin^6i7k) 



2 2 
"k-'^k 



2ek 



and 7k = ^(cosfc^; + cosfcj,). The magnon energy in the 
harmonic (or linear spin- wave) approximation is given by 
Ek- The field-induced planar anisotropy opens a gap at 
k = 0, while an acoustic branch is preserved for k — > Q. 

The lowest-order quantum correction to the magnon 
energy is determined by two types of anharmonicitics, 
cubic and quartic. The cubic term, which originates from 
the coupling between transverse and longitudinal fluctu- 
ations in the canted spin structure, is given by 



^3 = 



J sin 29 J2 e''*""' (4 + a.Oa]aj 



(7) 



/ H \ 

cos9[2Jsm9 - —j 



with rii — a|ai. The quartic term is 
Hi = j'S^i -- cos^9 [{ui + nj)aiaj +h.c.] 



^ L4 



(8) 



— niUj cos 29 — — sm^9 [al{ni + nj)aj + h. c] | . 

The frequency- independent l/S correction to the 
magnon energy due to Ji^ is most easily found by per- 
forming the Hartree-Fock decoupling in Eq. (|S]) with the 
following mean-field averages^^ 



4aj> = Xl'^kWk: (9) 



k k 

S = (flj ) = ^ Uki'k , A = (aiaj) = ^ 7kMkt'k ■ 

k k 

The frequency-dependent corrections from quartic anhar- 
monicitics appear only in the next order in l/S and are 
not considered here. 
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The role of the cubic anharmonicity is twofold. First, 
the Hartree-Fock decoupling in Eq. ([7]) yields quantum 
correction to the canting angle)^ 



sni a = sm f 



1 . 

l + -^(n-m- A) 



(10) 



Substituting sin^ back into Eq. (|4]) one obtains another 
frequency-independent correction to the spectrum. The 
combined Hartree-Fock correction, which includes the 
above two contributions, reads as 



A (1) 



A - n + sin26'(A + m)(l - 27k cos^O) (11) 
7k { n - A - m + cos 26* ( A sm^O + m cos^B) } 



with Wk being the dimcnsionless magnon energy, see 
Eq. ©. 

The second contribution is due to the remaining fluctu- 
ating terms in Tij, described by three-magnon processes: 



2! 



k,q 



^3 = ^E'^2(k,q) 



^Lq+Q^q^k + h. C. 



^k&H-k-q + h-C. 



k,q 



Decay and source vertices are given explicitly by 



(12) 
(13) 

(14) 



$i,2(k,q) = -i/cos0$i^2(k,q)/V25iV , 

$l(k,q) = 7k(Mk + Wk)(Wql'k-q+Q+ ^'q-Uk-q+q) 

+ 7q("q + «q)("kMk-q+Q+ I'kWk-q+q) 

+ 7k-q+Q(uk-q-|-Q+Uk-q+Q)(ukUq + l'kfq), 

$2(k,q) = 7k(Mk + Wk)(Wq?^k+q-Q+ WqUk+q-q) 

+ 7q("q + «q)(^tkWk+q-Q+ WkUk+q-q) 

+ 7k-|-q-Q(wk+q-Q+Uk+q-Q)(MkWq + UkUq). 

Next, we define the bare magnon Green's function 

Goi(k,a;) =a;-ek + jO (15) 

and calculate the second-order self-energy correction gen- 
erated by the cubic terms, which are represented by two 
diagrams in Fig. [2j 



Si(k,w) 
S2(k,c^) 



l*l(k,q) 



W - Cq - fk-q+Q + iO 



l*2(k,q)|2 



W -I- eq + Ek+q-Q - «0 



(16) 
(17) 



In the leading l/S order quantum correction is given by 
the on-shcU self-energy: 



Se-^^ = i;(k,ek) = -AJsin^e cos 



n 



l^l(k,q) 



2« (18) 
l^2(k,q)P 



^q+ ^k-q+Q— "^k Wk+ l^q+ Wk+q-Q 
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FIG. 3: (Color online) Magnon energy of the spin-1/2 square- 
lattice antiferromagnet below the decay threshold. Dashed 
lines show the harmonic spin- wave energy ek, solid lines repre- 
sent the spectrum with the lowest-order quantum correction. 
The inset in (a) defines the chosen path in the Brillouin zone. 



with the net renormalization for the excitation energy: 

ek = ek + Se^^^ + Se^^^ . (19) 
The magnon decay rate in the l/S order is given by 

Tk = ? E \^^(^' q)l''5(ek - eq - ek-q+q) . (20) 



Note, that the above Born approximation yields a spin- 
independent value of the decay rate Fk = 0(J). 

Although, our focus in Sec. IV is on the large-S* model, 
we now present numerical data for the spectrum renor- 
malization ([T^ in the S = 1/2 case. Similar results 
for the other values of spin can be obtained by rescaling 
quantum corrections in Figs.[3]and|l]by the factor 1/(25). 
For moderate magnetic fields, H < O.GHg, the magnitude 
of the correction to the magnon energy is comparable to 
that in zero fieldJ^i^ Numerical results for Ck along the 
symmetry directions in the Brillouin zone are shown in 
Fig. [3] for two values of external field. Contrary to the 
H = case, the renormalization is momentum-dependent 
with increasing deviations from the harmonic theory as 
the field approaches the decay threshold H* w 0.75-ffs- 
Specifically, there is a sizable dispersion along the mag- 
netic zone boundary X'X, which further increases in 
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FIG. 4: (Color online) Magnon energy and decay rate of 
the spin-1/2 square-lattice antiferromagnet above the decay 
threshold. Upper panels: dashed lines show the harmonic 
spin- wave energy £k, solid lines represent the spectrum with 
the lowest-order quantum correction tk- Lower panels: the 
decay rate Fk in the Born approximation. 



fields above one-half of the saturation field. This is pre- 
cisely the field regime where the hybridization between 
one- and two-magnon states grows substantially. 

The lowest-order spin- wave correction changes dramat- 
ically at higher fields. The renormalized magnon energy 
eic exhibits peculiar singularities in the form of jumps and 
spikes, see Fig. These anomalies signify a breakdown 
of the perturbative spin- wave expansion since the 1/5 
correction appears to be divergent for certain momenta. 
Such an outcome is rather surprising because the expan- 
sion parameter {a\ai) / {2S) decreases monotonically with 
the field and vanishes aX H = Hs- In particular, the 1/5 
spin-wave series for the SAFM converges rapidly for the 
ground-state energy and uniform magnetization^i^ which 
agree very well with the numerical results*^ 

The observed anomalies in the dynamical properties 
are clearly related to zero-temperature magnon decays: 
every singularity in ek is accompanied by a jump or by 
a peak in the decay rate T\^ }^-^^ The typical magnitude 
of Fk for high-energy magnons is comparable to or even 
exceeds ek, which means that the dynamical response in 
certain parts of the Brillouin zone is determined predom- 
inantly by the two-magnon continuum. A detailed con- 
sideration of the origin of the decay anomalies is provided 
in the next Section, while the asymptotic expressions for 



Fk-s-Q are derived in Appendix A. 



III. KINEMATICS OF MAGNON DECAYS 

The problem of spontaneous (zero-temperature) quasi- 
particle decays generated by cubic anharmonicities has 
a rather long history. The phenomenon is compre- 
hensively documented for phonons in crystals^ and 
for excitations in superfluid ^He<^i2^ Similar effects 
in quantum magnets have started to attract attention 
only recentlyj^ i^^'^^i"^^" — Appearance of spontaneous two- 
magnon decays for the SAFM is controlled by the energy 
conservation law: 



Ek — Eq + ek-q+Q 



(21) 



The ordering wave-vector Q enters the momentum con- 
servation condition because of the staggered canting of 
local moments If the energy conservation is satis- 
fied only for a trivial solution q = Q, then magnon with 
the momentum k remains stable. The decay rate Fk be- 
comes finite if there are non-trivial solutions of Eq. (|2ip . 
A less obvious but not less generic outcome of the mix- 
ing between one- and two-particle states is the trans- 
fer of some of the Van Hove singularities from the two- 
particle density of states onto the single-particle spec- 
trum producing various non-analyticities in the latter 
Thus, studying the kinematics of two-magnon decays for 
a given dispersion ek one should consider two related 
problems: (i) what is the decay region in the momen- 
tum space where excitations are unstable, and (ii) where 
do the renormalization corrections exhibit singular be- 
havior. In the present case, we are also interested in the 
evolution of both the decay region and the singularities 
as a function of external magnetic field. 

Let us begin with the decay region. At the boundary 
of such a region the single-magnon branch crosses with 
the bottom of the two-magnon continuum. Therefore, 
the boundary can be determined by solving the system 
of two equations: Eq. (j2ip and the extremum condition 
imposed on its r.h.s.: 



Vk-q+Q 



(22) 



where Vk = Vkek is a magnon velocity. Generally, 
one can envisage several types of solutions for the de- 
cay threshold i^i^ 

(i) decay threshold with the emission of an acoustic 
magnon q — > Q, determined by the condition 



Vk 



|vq| 



(23) 



(ii) decay threshold with the emission of two magnons 
with equal momenta q = i(k + Q) found by solving 



ek = 2e 



(k+Q)/2 



(24) 



and (iii) decay threshold with the emission of two 
nonequivalent magnons found from a direct numerical 
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FIG. 5: (Color online) Field evolution of the decay region of 
the SAFM in the Born approximation. Numbers correspond 
to values of H / Hs. 

solution of Eqs. ([?!]) and (P^ . Then, the decay region is 
the union of regions obtained for all three decay channels. 

Magnons in the square-lattice antiferromagnet remain 
stable up to the threshold field H* sa 0.7577^^^ Numeri- 
cal results for the decay regions in higher magnetic fields 
are shown in Fig.[S] In fields slightly above if*, the insta- 
bility first affects the acoustic branch producing a cigar- 
shaped decay region oriented along the TM diagonal in 
the Brillouin zone. Expanding energy of the acoustic 
magnon we obtain 

eQ+k«cfc + afc^, c = 2 V2 J5 cos e , 

a = -[t&v?e ' (25) 

where c is the acoustic magnon velocity and a is the coef- 
ficient of the leading nonlinearity in the dispersion, which 
depends on the momentum direction via the azimuthal 
angle if. Then, using smO = H/Hs, one can easily find 
that a changes its sign from negative in low fields to pos- 
itive in high fields. Such a change in convexity of the 
acoustic branch is generally known to result in the kine- 
matic instability with respect to two-particle decaysj^i^ 
see also discussion in Appendix A. The sign of a changes 
first for (y9 = 7r/4 at the threshold field 

H* = ^ Hs w 0.7559 Hs . (26) 
V 7 

The above consideration is based on the harmonic ex- 
pression for ek. Quantum corrections to the magnon 
dispersion may rcnormalizc the value p6|) . although the 
corresponding effect is expected to be small for S > 1/2. 

The described change in the curvature of the Gold- 
stone branch with increasing field is a rather general 
phenomenon. The acoustic mode changes convexity at 
the same field (pS)) in the cubic-lattice Heiscnberg an- 
tiferromagnet as well as in the stacked square-lattice 



model for any value of the antiferromagnetic coupling 
between layers. Moreover, the spin-wave velocity van- 
ishes at H = Hs for all Hcisenberg antiferromagnets. In 
this critical field the magnon branch is parabolic at low 
energies: eq+k oc k^. By the continuity argument, the 
spectrum preserves its positive curvature for a certain 
range of magnetic fields in the ordered phase, where the 
asymptotic form (j25p holds. This, in turn, implies that 
the two-particle decays are allowed at least for k — )• Q. 
A well-known example of such a behavior is provided by 
the hard-core Bose gas.— More generally, low-energy ex- 
citations in the ordered phase near XY quantum critical 
point are always unstable with respect to two-particle de- 
cays. Their decay rate has a universal form Fk oc k^^^^, 
in agreement with previous results for D = 2 and 3 .^^'^^ 
A more detailed discussion of the decay rate of low-energy 
magnons in the SAFM is provided in Appendix A. 

Magnetic systems with the 0(2) rotational symmetry 
about the field direction belong to the same XY univer- 
sality class as the Heisenberg SAFM in a field. An analy- 
sis of the kinematic conditions for magnon decays in the 
XXZ square-lattice antiferromagnet in a field parallel 
to the z-axis is briefly summarized in Appendix B. An- 
other common example is given by quantum spin systems 
with the singlet ground state and gapped triplet excita- 
tions. The field- induced ordering in such systems has 
been intensively studied both experimentally^^"— and 
theoreticallyi^"— The mechanism of spontaneous magnon 
decays discussed here applies equally to these magnets in 
the vicinity of the triplet condensation field He in the or- 
dered phase because of the duality between He and Hg 

For the Heisenberg SAFM, the decay region quickly 
spreads out across the Brillouin zone as the applied field 
increases above H* and at H > Q.9Hs most of the 
magnons become unstable already in the Born approx- 
imation, see Fig. [5] Up to ii « 0.85Hs the boundary of 
the decay region is entirely determined by the decay into 
a pair of magnons with equal momenta (|24|) . At higher 
fields H > 0.85Hs the decay channel with emission of an 
acoustic magnon (|23|) starts to prevail in some parts of 
the Brillouin zone producing a more complicated shape 
of the decay region, see Fig. [5l 

We now turn to the anomalous features in the magnon 
spectrum. A close inspection of Fig. [3] reveals two distinct 
types of singularities in it. The first type is characterized 
by a dip in ek and by a jump in Fk, while the anomaly 
of the second type consists in a jump in ek accompanied 
by a peak in Fk. The imaginary part of the magnon 
self-energy is related to the two-magnon density of states 
via Eq. (PP]) . Hence, in 2D the decay rate Fk exhibits a 
finite jump upon entering the continuum in accordance 
with the corresponding Van Hove singularity in the two- 
particle density of statesi^^ 

Fk ~ e(Afc) , Re Ek(ek) - In ^ , (27) 

\Ak\ 

where A is a cutoff parameter. The associated singularity 
in RcEk(ek) follows from the Kramcrs-Kronig relation. 
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less, the kinematic analysis above still predicts the cor- 
rect threshold field for magnon decays. This is because 
all higher-order n-particle decays with n > 3 are energet- 
ically forbidden if no two-magnon decays are allowed in 
the whole Brillouin mwe^ In our case, this means that 
magnons in the SAFM remain completely stable up to 
the two-magnon threshold field H* . Once two-magnon 
decays become possible at iJ > H* , the n-magnon de- 
cays are also kinematically allowed. Generally, we expect 
the corresponding decay rates to be significantly smaller 
than for the two-magnon decays because of their higher 
1/S' order. Moreover, the multi-magnon decays should 
not produce any singularities of the type ([27)) and (|28)) 
due to higher-dimensional integration in corresponding 
analogs of Eq. (|20l). 



FIG. 6: (Color online) Decay region and associated singulari- 
ties for H = 0.9-ffs. Line's type represents the decay channel; 
solid is for a pair of magnons with equal momenta, dashed is 
for a pair of magnons with different momenta, and dotted is 
for emission of an acoustic magnon. Line's color indicates the 
type of singularity in the continuum: blue is for the minimum 
and red is for the saddle-point. 

The second type of anomaly corresponds to an inter- 
section of the one-magnon branch with the saddle-point 
Van Hove singularity inside the continuum)^ 

Tk - In ^ , Re Sk(ek) - sign(Afc) . (28) 
|Afc| 

The boundary line given by Eq. (1^^ always corresponds 
to a local extremum and when it enters the interior of 
the decay region its type changes from a local minimum 
to a saddle point. The details of such a behavior are il- 
lustrated in Fig. ini Here, solid lines denote singularity 
contours for the decay into equivalent magnons, dashed 
lines indicate the decay threshold into a pair of magnons 
with different momenta, and dotted lines represent the 
decay threshold for emission of an acoustic magnon. The 
colors are used to distinguish whether the extremum cor- 
responds to a local minimum (blue) or to a saddle-point 
(red) in the continuum. The same singularity contour 
may correspond to a minimum in one part of the Brillouin 
zone and to a saddle-point in the other. This determines, 
in turn, whether the logarithmic singularity occurs in the 
real or imaginary part of the self-energy. Further en- 
hancement of the magnon decay rate may occur in the 
vicinity of the intersection of two saddle-point lines, an 
example which can be found in Fig. [B] along the TX line. 

We finish this section with a brief discussion of higher- 
order magnon decays. Similar to the consideration of 
the cubic nonlincarities in Sec. II, the quartic terms 
in Eq. ([5]), treated beyond the Hartree-Fock approxi- 
mation, produce decays in the three-particle channel: 
^q^P^k-q-p^k- In addition, higher-order bosonic terms 
omitted in Sec. II should open decays of one magnon 
into an arbitrary large number of magnons. Neverthe- 



IV. SELF-CONSISTENT THEORY OF 
MAGNON DECAYS 

A. One-magnon Green's function 

Singular quantum corrections to the magnon spectrum 
obtained in the leading 1/S' order are generally regular- 
ized by higher-order processes. Details of such a regular- 
ization depend on the specific shape of ek- An important 
aspect of the kinematics of the field-induced decays in 
SAFM is that the low-energy magnons created in a de- 
cay process are unstable by themselves. Qualitatively, 
the inverse life-time Fq of the decay products cuts off the 
singularities in (P7)) and ([^5]) by changing 

|Afc| ^ y^(Afc)2 + (Fq/Cq)2 , 

where Cq can be expressed through the velocities of the 
decay products. Accordingly, the self-consistent proce- 
dure, which replaces the bare magnon Green's functions 
in the self-energy diagrams in Fig. [2] with the renormal- 
ized ones, should produce non-singular quantum correc- 
tions. The main problem with such a self-consistent Born 
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FIG. 7: (Color onUne) Magnon decay rate Tk at H ^ 0.9Hs 
calculated in the Born approximation (dashed line) and in the 
iSCBA scheme for S = 1 (thin solid line, red) and for S — 5/2 
(bold solid line, blue). Dots in TM panel are — ImEk(ek) for 
5=1/2 from Ref. El. 



0.8H, 



II 

Co 



X' M X' M X' M X' M 




X' M X' M X' 



M X' 



m 



1 




rj 





0.14 




Co 



FIG. 8: (Color online) Intensity maps of the magnon decay rate Fk calculated within the iSCBA for several values of external 
magnetic field and for the two values of the spin: 5=1 (upper panels) and 5 = 5/2 (lower panels). 



approximation (SCBA) is that it opens an unphysical 
gap for acoustic magnons in violation of the Goldstone 
theorem. Here we avoid this problem by performing a 
restricted self-consistent calculation, which takes into ac- 
count the imaginary part of the magnon self-energy and 
neglects the real part of it. Such an approximation, which 
is referred to as iSCBA in the following, is expected to 
yield reasonable accuracy for spins 5 > 1. Indeed, in the 
low-field region the nonsingular quantum corrections to 
the magnon energy do not exceed 15-18% in the case of 
S = 1/2, see Sec. II. For larger spins this correction is 
further reduced by a factor 1/(25') and already for 5 > 1 
it corresponds only to a small shift of the magnon energy. 

In the simplest realization of the iSCBA we further ne- 
glect the frequency dependence of ImSk(a;) and impose 
the following form of the one-magnon Green's function: 



G ^ (k, cj) = cj — Ek + iFi, 



(29) 



The magnon decay rate Ft is calculated self-consistently 
from 



Tm S](k, Ek) 



(30) 



with S(k, Ek) from Fig. [2l Physically, the approximation 
amounts to assuming the Lorentzian shape of the 
quasiparticle peak in the dynamical response. Such a 
self-consistent scheme has been previously applied to the 
problem of phonon broadening in superfluid helium^ and 
in quasicrystals42 In the present case, with two types of 
cubic vertices the self-consistent equation on Fk becomes 



Fl 



2 ^ 
q 



|a>l(k,q)|2(Fq + Fk-q+Q) 



(Ck - Eq - Ck-q+q)^ + (Fq + Fk-q+Q)^ 
|$2(k,q)|2(Fq + Fk+q-Q) 



(ek + Cq + Ek+q-q)^ + 



Tk+q-q)^ 



(31) 



One important consequence of the self-consistent con- 
sideration is that spontaneous decays are not restricted 
to the decay region in Fig. [5l defined in the Born approxi- 
mation. Since the finite decay rate of magnons inside the 
decay region means an uncertainty in their energy, the en- 
ergy conservation condition (|2ip is relaxed and magnons 
just outside the decay region are also allowed to decay. 
Therefore, the decay threshold boundary is smeared and 
the decay probability disappears gradually rather than 
in a step-like fashion. 

Numerical solutions of Eq. ((3T|) have been obtained 
iteratively using two different grids in the momentum 
space: a sparse mesh for the momentum k of 40 x 40 
points in the full Brillouin zone and a much tighter 
grid for integration over the momentum q with up to 
1000 X 1000 points. The values of Fq for the integrand 
were obtained from Fk by a bicubic interpolation. To 
improve convergence at low energies we also impose the 
asymptotic form Fq+k ^ fc"^ for small fc, see Appendix A. 
Sufficient accuracy is typically achieved after 5-8 itera- 
tion steps. Note, that the second term on the r.h.s. of 
Eq. ([5T|) from the source diagram in Fig.[5]is always small 
and can be safely neglected. 

Figure [7] compares Fk obtained in the Born approxi- 
mation with the iSCBA results for 5 = 1 and S = h/2 
&i H = 0.9Hs. One can see that the singular behav- 
ior of Fk near the Van Hove singularities is completely 
removed in the self-consistent calculation. The typical 
amplitude of the decay rate in Fig. [7] away from F and 
M points is Fk ^ 0.5 J, which should be compared with 
the typical magnon energies ek ^ 4:JS. Another notable 
aspect of the self-consistent regularization is in the spin- 
dependence of the decay rate: Fk is larger for larger spins 
in most of the Brillouin zone. Such a dependence should 
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not be confused with the relative strength of the spec- 
trum broadenmg Fk/ek ^ Tk/JS*, which is stih larger 
for smaher values of the spin. The decay rate approaches 
the spin-independent Born result only in the limit 

5 — >■ oo. In particular, while Fk — rj^"^ + 0(1/ S) away 
from the singular points, one finds Fk ~ In 5 in their 
vicinity 1^ 

Further details on the momentum distribution of the 
magnon decay rate in different fields arc provided in the 
form of intensity maps for Fk/4JS' in Fig. |51 Compari- 
son of the upper {S = 1) and the lower (S = 5/2) row 
reveals similar patterns in the momentum distribution 
of Fk for these two values of spin. The intensity of the 
decays in Fig. |8] is determined by cither the volume of 
the decays or by the proximity to the Van Hove singu- 
larities of the two-magnon continuum. Specifically, in 
the lower fields H < 0.9 Hs the magnon damping is most 
significant in the broad region around (7r/2,7r/2) point. 
This is a consequence of the large phase space volume 
for the two-particle decays in this region. In higher fields 
the maximums in the decay rate shift towards the TX 
and rX' lines. The enhancement of Fk here is due to 
an intersection of the two singularity lines discussed in 
Sec. III. Overall, the maximum magnon decay rate in all 
fields docs not exceed T^^^ ^ 0.7-0.8J. Thus, we may 
conclude that the spin waves in large-S' quantum anti- 
ferromagnets above the decay threshold H* are damped 
but still well-defined quasiparticles. 

The above conclusion seems to be at odds with the be- 
havior of the spin- 1/2 SAFM. In this case, magnons were 
shown to be overdamped in most of the Brillouin zone 
by both the spin-wave calculation, which used a differ- 
ent version of the SCBA^ and by the numerical QMC 
simulations Such a difference in the effect of magnon 
interaction between 5 = 1/2 and S > 1 models deserves 
a separate discussion. We compare directly the imagi- 
nary part of the magnon self-energy at the quasiparti- 
cle pole, — Iml]k(ek), for the 5=1/2 case obtained in 
Ref. [l^ with Fk found in this work, see the data points 
along TM direction in Fig. [71 One can see that 5 = 1/2 
and 5 > 1 results are quantitatively close. Naively, this 
would mean that quasiparticles in the 5 = 1/2 antiferro- 
magnet should also be reasonably well-defined. However, 
this consideration neglects the spectral weight redistribu- 
tion. The analytical calculation^^ and the QMC studyii 
have shown a significant non-Lorentzian broadening of 
the spectral lines, which takes the form of a double-peak 
structure in the dynamical susceptibility. Therefore, the 
overdamping in the 5 = 1/2 case is produced by two 
effects: broadening of the quasiparticlc peak and spec- 
tral weight redistribution. Note, that the latter effect 
was completely excluded in the recent hard-core boson 
studyr^ which explicitly neglected the frequency depen- 
dence of the magnon self-energy. Such an approximation 
may have lead Ref. [l^ to the conclusion of only weak- 
to-moderate damping of spin excitation in the 5 = 1/2 
SAFM in a field, in contrast with Refs. [H andllTl. 

Since the iSCBA developed here also neglects the spec- 



tral weight redistribution, one can question its validity 
on the same grounds. To verify the validity of this ap- 
proach we have checked the accuracy of the Lorentzian 
approximation used in (j29p by keeping the full frequency 
dependence of Iml]k(w) and by utilizing general repre- 
sentation for the magnon Green's function 



G(k,w) 



da; ■ 



Ak{x) 



(32) 



The self-consistent equation is formulated for the spec- 
tral function j4k(w) excluding again the real part of the 
magnon self-energy. Numerical results (not shown) ex- 
hibit only a small asymmetry of magnon peaks and only 
for the lowest 5 = 1 case. This justifies the use of the 
frequency-independent iSCBA scheme for 5 > 1 antifer- 
romagnets and confirms the perturbative role of damping 
for them once the singularities are regularized by a self- 
consistent calculation. 



B. Dynamical structure factor 

The inelastic neutron scattering experiments probe the 
spin-spin correlation function 



dt 
2^ 



(5£(t)5: 



(33) 



defined in the global laboratory frame, while the spin- 
wave calculations are conveniently performed in the local 
rotating frame, see Sec. HIl First, we relate these two 
forms of the correlation function with the help of Eq. ^ : 

5^°-^°(k,w) = sin26i5^^(k,w)+cos26'5""(k-Q,a;), 
5^°^°(k,w) = cos^ 61 5"^^ (k -Q,w) + sin2 6* 5"^ (k,cj), 
5^°«°(k,w) = 5^^(k,cj) . (34) 

Here we have omitted the cross-terms [a ^ /3) on the 
r.h.s., which can be shown to be small numerically. 
Second, we relate the dynamical structure factor (p4)) 
to the time-ordered spin Green's function Q"'^ (k, t) = 
—i{TS^{t)S^-^ using the fluctuation-dissipation theo- 
rem at zero temperature: 



5"'3(k,cj) = --Img"'^(k,a; 



1 



(35) 



Finally, we express spin operators via the Holstein- 
Primakoff bosons and, after some of algebra, obtain the 
following expressions 

5^^(k,w) = -5A+ (wk + wk)^ ImG'(k,w) , 

5^^(k,cj) = -5A_ (uk-wk)^ ImG(k,w) , (36) 

5"'^(k,w) = --^^ (itk-qi'q + UqUk-q)^ 

q 

X Im < / da; G(q, a;) G(k - q, w - x) > . 
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FIG. 9: (Color online) Transverse part of the dynamical struc- 
ture factor (gray scale) computed using the iSCBA for 5 = 1 
and for the two values of external magnetic field. Dashed 
lines (blue) enclose the region ek ± Tk corresponding to the 
full- width at half maximum of the corresponding quasiparticle 
peaks. 



Here A± = 1 — (2n ±5)123 are the Hartree-Fock spin re- 
duction factors and G(k, uj) is the magnon Green's func- 
tion from ((29)) . 

In the equation above, spin fluctuations are separated 
into transverse {xx^ yy) and longitudinal {zz) compo- 
nents relative to the local spin direction, which corre- 
spond to scattering from states with odd and even num- 
ber of magnons. respectively An important effect of 
the spin canting on the dynamical structure factor is the 
redistribution of the spectral weight over the two well- 
separated transverse modes^^ which overlap only at the 
magnetic zone boundary. As a result, Eq. ([34]) contains 
a mixture of two transverse contributions at momenta k 
and k — Q. The first mode corresponds to spin fluctu- 
ations perpendicular to the magnetic field direction (in- 
plane mode), while the second mode, shifted by the order- 
ing wave- vector Q, corresponds to the fiuctuations along 
the field direction (out-of-plane mode). In strong mag- 
netic fields the in-plane mode is enhanced and dominates 
over the out-of-planc fluctuations, which gradually disap- 
pear as H ^ Hs- Note also the existence of two distinct 
longitudinal multimagnon continua associated with each 
of the transverse branches. 

Results summarized in Eqs. ([M]) and ([55]) are quite 
general and can be adapted to a variety of specific ex- 
perimental situations. As a particular example, the dy- 
namical structure factor accessible to unpolarized neu- 
tron scattering experiment with a vertical magnetic field 
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FIG. 10: (Color online) Same as Fig.[9]for S = 5/2. 



reads as 5(k,w) = 5_L(k,w) + S'"""'{k,uj). The explicit 
expression for 5^(k, oj) depends on the experimental de- 
tails via a momentum-dependent polarization factor. In 
the following, such a factor is deliberately ignored and 
the expression 5(k, w) = iS^"^" + 5^0^° + S''°''° is used 
for an illustration. The combined 5(k, a;) computed us- 
ing the self-consistent magnon Green's function obtained 
in previous subsection is shown for S = 1 and S* = 5/2 
antiferromagnets in Figs. and [TUl respectively. One can 
clearly see the presence of two transverse modes with the 
more intense in-plane and the weaker out-of-plane mode. 

A distinct fingerprint of spontaneous magnon decays 
as compared to other possible damping mechanisms is 
their characteristic dependence on the momentum and 
on the applied magnetic field. In the field range H* < 
H < 0.9Hs the magnon damping is most significant along 
the magnetic zone boundary X'X with the maximum at 
(7r/2,7r/2). In higher magnetic fields H > 0.9Hs the 
strongest damping occurs along the TX line, where the 
singularity-enhanced decays are most pronounced. Our 
predictions for the former field regime are in a quali- 
tative agreement with the very recent inelastic neutron 
scattering results of Masuda et for the spin-5/2 lay- 
ered square-lattice antiferromagnet Ba2MnGe207. How- 
ever, the largest magnon linewidth measured experimen- 
tally exceeds our theoretical estimate 2rk — 1.5J for 
H = 0.82Hs by a factor of 2-3. The larger experi- 
mental linewidth may be related to the overlap of the 
two transverse modes along the magnetic zone bound- 
ary. Therefore, we propose that the inelastic neutron 
scattering measurements at H > 0.9Hs in the vicinity of 
(7r/2,0) point where the two transverse modes are well 
separated may provide a benchmark for the field-induced 
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spontaneous magnon decays in the SAFM. 



V. CONCLUSIONS 

We have presented a detailed analysis of the field- 
induced decay dynamics in the square-lattice antifer- 
romagnct within the framework of the spin-wave the- 
ory. In magnetic fields exceeding the threshold field for 
two-particle decays, the 1/5 quantum correction to the 
magnon spectrum exhibits singularities for certain mo- 
menta. These are related either to the decay thresholds'^ 
or to the saddle-point Van Hove singularities in the two- 
magnon continuumj^ Such singularities invalidate the 
usual 1 / 5'-expansion because the resultant rcnormalizcd 
spectrum contains divergences. The encountered prob- 
lem is rather generic and must be common to a variety 
of 2D models regardless of the on-site spin value. For 
the systems with large spins the situation is especially 
aggravating as the 1/S expansion turns from a reliable 
approach to the one producing unphysical divergences, 
which complicate any sensible comparison with experi- 
mental or numerical data. 

In the present work wc have developed a self-consistent 
regularization scheme, which is applicable to a variety 
of problems in the decay dynamics of large-S" quantum 
antiferromagnets. The self-consistently calculated damp- 
ing in the SAFM is free from singularities and gives an 
upper limit on the decay rate: (Fk/ek) ^ 0.2/5. Over- 
all, spin-waves in the SAFM with S* > 1 in magnetic 
field H > H* do acquire finite broadening, but remain 
well-defined quasiparticles. The regularized singularities 
also lead to a parametric enhancement of the magnon 
damping and produce a distinct field- and momentum- 
dependence of Fk, which is illustrated for the 5 = 1 and 
5 = 5/2 in Figs. [TlfTUl These results can be used for 
quantitative comparison in future experimental studies 
of the field-induced spontaneous magnon decays. 

For the spin- 1/2 SAFM the situation appears to be 
somewhat more delicate. Previous analyticali^ and 
numericalii studies have predicted overdampcd one- 
magnon excitations in a large part of the Brillouin zone. 
The SCBA scheme used in Ref. [l^ includes a self- 
consistent renormalization of only one inner magnon line 
in the decay diagram of Fig. [2] and is, in a sense, not 
as consistent as the present approach. On the other 
hand, that approach does take into account the real part 
of the spectrum renormalization and, as a consequence, 
the quasiparticle weight redistribution. Such an effect 
can be deemed small for larger spins, but it is more im- 
portant for 5 = 1/2 and is likely to contribute to fur- 
ther enhancement of the damping in this case. Since 
the imaginary part of the self-energy at the quasiparticle 
pole for 5 = 1/2 in Ref. [l^ correlates closely with the 
iSCBA results of the present work (Fig. [7]), we believe 
that this is the correct explanation of the differences be- 
tween 5=1/2 and 5 > 1 cases. There are also qualita- 
tive and quantitative similarities of the analyticali^ and 



numerical resultSfii in particular the double-peak struc- 
ture in the spectral function. While the QMC approach 
may involve its own uncertainties due to numerical in- 
terpolation from imaginary to real frequencies such 
a correspondence between results of two very different 
methods is encouraging. Further theoretical efforts may 
be needed to clarify completely the detailed behavior of 
the dynamical structure factor for the spin-1/2 SAFM. 
Finally, performing inelastic neutron scattering measure- 
ments on suitable spin-1/2 compounds would also be im- 
portant. 
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Appendix A: Decay of low-energy magnons 

In the case of the Heisenberg SAFM and other two- 
sublattice antiferromagnets in applied magnetic field the 
low-energy spectrum consists of the single weakly nonlin- 
ear acoustic branch In the following consideration 
all momenta arc taken relative to the magnetic ordering 
wave- vector Q: eq+k — > fk ~ cfc + ak^ . First, wc verify 
explicitly the energy conservation condition for the two- 
particle decays (|2ip . If the nonlinearity of the spectrum 
is weak ak^ ^ cfc, the two magnons emitted in a sponta- 
neous decay process have their momenta q and q' = k— q 
almost parallel to the direction of k. Then, 



|k-q| « fc-g- 



kqip'^ 



2{k-q) ' 



(Al) 



where is a small angle between q and k, and the energy 
conservation can be written as 



ckq 



3akq{k — q) 



ckqtf^ 
2(fc - q) 



(A2) 



2{k-q) 

One can sec that the nontrivial solution Lp ^ Lp^ exists 
only for the positive sign of the cubic nonlinearity 

Next, we derive the asymptotic, k,q ^ form for the 
decay vertex $i(k, q), see Eq. (|14p . In order to obtain 
the correct expression one needs to retain three leading 
terms in the expansion of the Bogolyubov coefficients Uk 
and Vk in small k. This yields 



$1 



cos 9 



1/2 



i/fcg(fc - q) 



\\J kq{k - q) 



4 cos^ 6 



c 



(A3) 
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For momenta k and q, which satisfy the energy conser- 
vation condition (pij) . the decay vertex further simpUfies 
into 



1.1 (k,q) 



cos 6 



4 cos2 9 \ y/2 



1/2 



Vkq{k - q) . (A4) 



Substituting (jA2p and (jA4p into the Born expression for 
the decay rate in Eq. pO|) one obtains 



3J 
16^ 



tan^6' 



(A5) 



The bare values of c and a are given by Eq. (^5)) . al- 
though one may generally use renormalized parameters. 
Note, that the asymptotic expression (|A5p differs from 
the result provided in Ref. Il9l . Therefore, we have ver- 
ified that our expression for Fk in (jASp agrees with the 
direct numerical integration of Eq. (|20l) in the limit of 
small k. 

As H Hs, the velocity of the acoustic mode de- 
creases and the condition of a weak nonlinearity applies 
to a progressively narrower range of momenta <C 
8 cos^ 9 reducing the range of validity of the asymptotic 
expression (|A5I) . Outside that domain one can use the 



parabolic form of the magnon dispersion ek ~ JSk to 
derive another useful asymptotic expression for the de- 
cay rate. In this regime Wk ~ 1 and Wk = 0{cos9) « 
and the decay vertex becomes $i(k, q) w —2. The 
angle between the emitted and the initial magnon can 
now be large: cost^o = l/k. Performing an analyti- 
cal integration in Eq. (PO)) and taking into account that 
cos^ 9 « 2(1 — H/Hs) we obtains 



16 1 



H 



(A6) 



which is valid for 4^1 - H/Hs < fc < 1. 

Modifications to the asymptotic expression (jASp are 
also expected in the field regime just above the decay 
threshold field, H — > H* . In this case, a — ?> and 
the coefficient in front of k^ in the Born expression for 
Fk in (|A5P diverges. Such a nonanalytic behavior is 
nothing but the long-wavelength version of the decay 
threshold singularities. The self-consistent regularization 
([5T|) of Sec. IIVI should remain applicable in this limit. 
In order to verify that the power-law behavior (jA5P is 
not modified within the SCBA we substitute a general 
power-law ansatz for the decay rate Fk = and as- 
sume that the damping is much larger than the nonlin- 
earity but still much smaller than the magnon energy: 
ak^ <^ f3k^^ ^ ck. In this case, the decay angle is scaled 
as (p^ ^ fc"^^ with q ^ k, which makes the decay vertex 
$i(k, q) ^ instead of (jA4p . The power counting 

on both sides of Eq. ([51]) yields a unique solution: n = 3. 
Therefore, the Born exponent Fk = Pk^ is not changed 
by the self-consistent procedure. The damping coefficient 
/3 does not diverge near the decay threshold field anymore 
because a drops out from the self-consistent equation on 



Co 



, 1 I 1 I 1 I 1 


, 1 < 1 < 1 < 1 


easy-plane 


^\ easy-axis 


















. 1 . 1 . 1 . 1 


. 1 . 1 . 1 . 1 



1 



0.4 0.6 

e 



0.6 0.4 

1/e 



FIG. 11: (Color online) Phase diagram of the XXZ model 
in applied field as a function of anisotropy parameter e. The 
field range where spontaneous magnon decays are allowed is 
shown in gray. 



Fk. Instead, /3 exhibits a step- like behavior in (H — H*), 
which is in an accord with the 2D character of the Van 
Hove singularity at the border of the two-particle contin- 
uum. 



Appendix B: Decays in XXZ SAFM 

The Hamiltonian of the XXZ antiferromagnet on a 
square-lattice is 

The values of the anisotropy parameter < e < 1 de- 
scribe the spin system with the easy-plane anisotropy, 
whereas e > 1 corresponds to the easy-axis case. 

Minimization of the classical energy yields the transi- 
tion field into a fully saturated state 



Hs = 4J5(1 + £) 



(B2) 



For the easy-axis case, the spin-fiop transition between 
the coUinear state and the canted antiferromagnetic state 
takes place at 



HsF = 2JS^2ie^ - 1) 



(B3) 



The harmonic spectrum in the canted antiferromag- 
netic state is given by the same expression Eq. ^ as 
for the Heisenberg SAFM with the substitution cos 29 — )■ 
(e cos^ 9 — sin^ 9). Performing the same type of kinematic 
analysis as in Sec. Ill, we find that the curvature of the 
acoustic branch changes at the threshold field 



H* = 4JS 



(B4) 



Above this field spontaneous two-magnon decays are 
kinematically allowed. The phase diagram in the H-e 
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plane is summarized in Fig. [TTJ Note, that in the right 
panel the units of the field are different from the left 
panel. 

Several observations are in order: (i) for arbitrary 
e there is a finite range of fields below the saturation 
field where magnons are unstable; (ii) the spin-flop fleld 
HsF in the easy-axis regime is always below the de- 
cay threshold field i/*; (iii) the easy-plane anisotropy 
pushes the decay instability to lower fields compared to 
the isotropic model. In the XY-limit (e = 0) the ratio 
H* /Hs = \l\fl 0.38 is two times smaller than for the 



Heisenberg antiferromagnet. This might be important 
for the search of experimental systems where the phe- 
nomenon of magnon decays can be observed. Note also 
that the easy-plane anisotropy appears naturally in the 
mapping of the BEC-type transition of the triplet excita- 
tions in quantum magnets with the singlet ground state 
onto the saturation field transition of an effective pseudo- 
spin-1/2 model The typical value of the anisotropy pa- 
rameter in these problems is e = 0.5, which implies that 
spontaneous magnon decays must exist in an extended 
interval of fields above the condensation field H^. 
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